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We studied the influence of the trigonal distortion of the regular octahedron along the (111) 
direction, found in the C0O2 layers. Under such a distortion the t^ g orbitals split into one a\ g and 
two degenerated e' g orbitals. We focused on the relative order of these orbitals. Using quantum 
chemical calculations of embedded clusters at different levels of theory, we analyzed the influence 
of the different effects not taken into account in the crystalline field theory ; that is metal-ligand 
hybridization, long-range crystalline field, screening effects and orbital relaxation. We found that 
none of them are responsible for the relative order of the ti g orbitals. In fact, the trigonal distortion 
allows a mixing of the ti g and e g orbitals of the metallic atom. This hybridization is at the origin 
of the ai g -e' g relative order and of the incorrect prediction of the crystalline field theory. 



I. INTRODUCTION 

Since the discovery of super-conductivity in the hy- 
drated Nao.3sCo02 — I.3H2O 1 compound and of the very 
large thermopower in the Nao.ri^CoO^ members of the 
same family, the interest of the community in systems 
built from C0O2 layers has exploded. The first step in 
the understanding of the electronic properties of transi- 
tion metal oxides, such as the Co02-based compounds, 
is the analysis of the crystalline field splitting of the d 
orbitals of the transition metal atom. Indeed, depending 
on this splitting, the spin state of the atom, the nature of 
the Fermi level orbitals, and thus the Fermi level proper- 
ties will differ. 

The C0O2 layers are built from edge-sharing C0O6 oc- 
tahedra (see figure [T]). In these layers, the first coordina- 




FIG. 1: Schematic representation of the C0O2 layers. 

tion shell of the metal atom differs from the regular octa- 
hedron by a trigonal distortion along the three- fold (111) 
axis (see figureE]). In all known materials (whether cobalt 
oxides or other metal oxides such as L1VO2, NaTiC>2, 
NaCr02, etc. . . ), this distortion is in fact a compression. 
The local symmetry group of the metal atom is lowered 
from Oh to £*3d- The Ti g irreducible representation of 



the Oh group is thus split into one E g and one A\ g rep- 
resentations. The relative energies of the resulting e' g 
and aig orbitals (see figure [6} has been a subject of con- 
troversy in the recent literature, as far as the low spin 
Co 4+ ion is concerned. At this point let us point out 
the crucial importance of the knowledge of this energetic 
order for the understanding of the low energy properties 
of the C0O2 layers. Indeed, the possible existence of an 
orbital order, as well as the minimal model pertinent for 
the description of these systems depend on this order. 

Authors such as Maekawa-, following the crystalline 
field theory, support that the a\ g orbital is of lower en- 
ergy than the two degenerated e g ones, leading to an 
orbital degeneracy for the Co 4+ ion. On the contrary, 
ab initio calculations, both using periodic density func- 
tional methods* and local quantum chemical methods 
for strongly correlated systems*, yield an a\ g orbital of 
higher energy than the e' g ones, and a non degenerated 
Fermi level of the Co 4+ ion. Angle Resolved Photoemis- 
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FIG. 2: Schematic representation of cobalt 3d splitting. 6 
represents the angle between the z axis — the 3-fold (111) 
axis of the CoOe octahedron — and the Co — O direction. 
80 = arccos(^^=*J ~ 54.74° is the 8 angle for the regular 
octahedron. 

sion Spectroscopy (ARPES) experiments were performed 
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on several C0O2 compounds 6 . This technique probes the 
Fermi surface and clearly shows that the Fermi surface 
of the C0O2 layers is issued from the a\ g orbitals, and 
not at all from the e' g orbitals (orbitals of E g symme- 
try, issued from the former t 2g orbitals), supporting the 
ab-initio results. 

In the present work, we will try to understand the rea- 
sons why the crystalline field model is unable to find the 
good energetic order of t 2g orbitals in such trigonal dis- 
tortions. Several hypotheses can be made to explain the 
orbital order : the derealization of the metal 3d orbitals 
toward the ligands, the fact that the electrostatic poten- 
tial of the whole crystal differs from the one assumed in 
the crystalline field model, the correlation effects within 
the 3d shell, the screening effects, etc. All these hypothe- 
ses will be specifically tested on the Co 4+ (3d 5 ) ion that is 
subject in this work to a more thorough study than other 
metal fillings. Nevertheless, other metal fillings (3d 1 to 
3d 3 , that can be found in vanadium, titanium chromium, 
. . . oxides) will also be studied. We will see the crucial 
importance of the band filling on the t 2g orbitals order. 
In this work we will focus only on the Oh to trigonal 
distortion, subject of the controversy. 

The next section will present the method used in this 
work, section three and four will reports the calculations 
and analyze them, finally the last section will be devoted 
to the conclusion. 



II. COMPUTATIONAL METHOD AND 
DETAILS 

The energy of the atomic 3d orbitals is an essentially 
local value, as supposed in the crystalline field model. 
However its analysis exhibits some non local contribu- 
tions. Indeed, orbitals energies can be seen as resulting 
from the following terms: 

• the electrostatic potential due to the first coordi- 
nation shell — in the present case, the six oxygen 
atoms of the octahedron, further referred as nearest 
neighbor oxygens (NNO) — , 

• the electrostatic potential due to the rest of the 
crystal, 

• the kinetic energy that includes the hybridization 
of the metal orbitals with nearest neighbor ligands, 

• the Coulomb and exchange contributions within 
the 3d shell, 

• the radial relaxation of the 3d orbitals, 

• and finally the virtual excitations from the other 
orbitals that are responsible for the screening ef- 
fects. 

All these contributions, excepts for the electrostatic po- 
tential due to the rest of the crystal (nucleus attractions 



and Coulomb interactions), are essentially local contri- 
butions^ and known to decrease very rapidly with the 
distance to the metal atom. In fact, they are mostly re- 
stricted to the first coordination shell of the cobalt. On 
the contrary, the Madelung potential retains the result- 
ing non local contributions from the nucleus attraction 
and the Coulomb electron-electron repulsion. It is known 
to be very slowly convergent with the distance. We thus 
made calculations at different levels, including first all 
the above effects, and then excluding them one at the 
time, in order to end up with the sole effects included in 
the crystalline field model. 

The calculations will thus be done on CoOg or Co 
fragments. Different embedding and different levels of 
calculation will be used. The Co — O distance will be 
fixed to the value of the super-conducing compound, i.e. 
R Co _o = 1.855 A. The angle 9 between the Co - O di- 
rection and the z axis (see figure 15]) will be varied from 
to 90°. 

The calculations will be done at the Complete Active 
Space Self Consistent Field + Difference Dedicated Con- 
figurations Interaction^ (CASSCF+DDCI, see subsec- 
tion [TTX| level for the most involved case, using the core 
pseudopotential and basis set of Barandiaran et alJ^. 
The fragment used will include all the first coordination 
oxygens in addition to the cobalt atom. The embed- 
ding will be designed so that to properly represent the 
full Madelung potential of the super-conducting mate- 
rial, and the exclusion effects of the rest of the crystal on 
the computed fragment electrons (see reference 5 - for fur- 
ther details). For the simplest case a minimal basis set 
derived from the preceeding one will be used and only 
the cobalt atom will be included in the computed frag- 
ment. The effect of the crystalline field will be described 
by —2 point charges located at the positions of the first 
coordination shell oxygens. The calculations will be done 
at the CASSCF level only. Between these two extreme 
cases, several intermediate ones will be considered, in or- 
der to check the previously enumerate points. 

The electrostatic potential due to the cobalt first oxy- 
gen neighbors (NNO) , as well as the unscreened Coulomb 
and exchange contributions within the 3d shell, are in- 
cluded in all calculations. The electrostatic potential is 
treated either through the inclusion of the NNO in the 
computed fragment or through —2 point charges. The 
Coulomb and exchange contributions are treated through 
the CASSCF calculation. The electrostatic contribution 
of the rest of the crystal is included only in the most 
involved calculations, using an appropriated embedding 
of point charges and Total Ions pseudo-Potentialii. The 
hybridization of the metal 3d orbitals is treated by in- 
cluding cxplicitcly the NNO in the considered fragment 
(CoOg). The radial relaxation of the 3d orbitals is treated 
when extended basis set are used. When a minimal basis 
set is used, the radial part of the orbitals is frozen as in 
the high spin state of the isolated Co 4+ ion. Finally, the 
screening effects are treated only when the calculation is 
performed at the CASSCF+DDCI level. 
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A. The CASSCF and DDCI methods 

Let us now described shortly the CASSCF and DDCI 
ab initio methods. These methods are configurations 
interaction (CI) methods, that is exact diagonaliza- 
tion methods within a selected set of Slater's deter- 
minants. These methods were specifically designed to 
treat strongly correlated systems, for which there is no 
qualitative single-determinant description. The CASSCF 
method treats exactly all correlation effects and exchange 
effects within a selected set of orbitals (here the 3c? shell of 
the cobalt atom). The DDCI method treats in addition 
the excitations responsible for the screening effects on 
the exchange, repulsion, hopping, etc. integrals. These 
methods are based on the partitioning of the fragment 
orbitals into three sets 

the occupied orbitals that are always doubly- 
occupied in all determinants of the Complete 
Active Space or CAS (here the cobalt inner 
electrons and the NNO ones), 

the active orbitals that can have all possible occupa- 
tions and spins in the CAS (here the cobalt 3d or- 
bitals), 

the virtual orbitals that are always empty in the 
CAS. 

The CASCI method is the exact diagonalization within 
the above defined Complete Active Space. The CASSCF 
method optimizes in addition the fragment orbitals in or- 
der to minimize the CASCI wave function energy. This 
is a mean-field method for the occupied orbitals but all 
the correlation effects within the active orbitals are taken 
into account. Finally the DDCI method uses a diago- 
nalization space that includes the CAS, all single- and 
double-excitations on all determinants of the CAS, ex- 
cept the ones that excite to occupied orbitals into two 
virtual orbitals. Indeed, such excitations can be shown 
not to contribute — at the second order of perturba- 
tion — to the energy differences between states that dif- 
fer essentially by their CAS wave function. Therefore, 
they have little importance for the present work. The 
DDCI method thus accurately treats both the correla- 
tion within the CAS and the screening effects. 

Compared to the very popular density functional meth- 
ods, the CAS+DDCI method presents the advantage of 
treating exactly the correlation effects within the 3d shell. 
This is an important point for strongly correlated mate- 
rials such as the present ones. Indeed, even if the DFT 
methods should be exact provided the knowledge of the 
correct exchange-correlation functional, the present func- 
tionals work very well for weakly correlated systems, but 
encounter more difficulties with strong correlation effects. 
For instance the LDA approximation finds most of the 
sodium cobaltites compounds ferromagnetic 4 in contra- 
diction with experimental results. LDA+U functionals 
try to correct these problems by using an ad hoc on-site 



repulsion, U, within the strongly correlated shells. This 
correction yields better results, however it treats the ef- 
fect of the repulsion within a mean field approximation, 
still lacking a proper treatment of the strong correlation. 
The drawbacks of the CAS+DDCI method compared to 
the DFT methods are its cost in term of CPU time and 
necessity to work on formally finite and relatively small 
systems. In the present case however, this drawback ap- 
pear to be an advantage since it decouples the local quan- 
tities under consideration from the dispersion problem. 



III. RESULTS AND ANALYSIS 

Let us first attract the attention of the reader on what 
is supposed to be the energy difference between the e' g 
and dig orbitals of the Co 4+ ion in an effective model. 
In fact, the pertinent parameters for an effective model 
should be such that one can reproduce by their means 
the exact energies or, in the present case, the ab-initio 
calculation of the different Co 4+ atomic states. It results, 
that within a Hubbard type model, the pertinent effective 
orbital energies should obey the following set of equations 

£(|a lg )) - As{e' g )+e(a lg ) + 2U + 8U' - U H 
E(\e' s )) = 3e(e' g ) + 2e(oi fl ) + 2U + W - U H 
AE = E(\e' s ))-E(\ ais }) 

= e(ais) - e ( e ' g ) 

where the schematic picture of the |e' ) and |ai g ) states 
is given in figure [3J £(e' g ) an d £(oig) are the effective 
orbital energies of the e' g and a\ g atomic orbitals, U is 
the effective electron-electron repulsion of two electrons 
in the same cobalt 3d orbital, U' the effective repulsion 
of two electrons in different cobalt 3d orbitals and Jh 
the atomic Hund's exchange effective integrals within the 
cobalt 3d shell. 
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FIG. 3: Schematic representation of the Co 4+ states of in- 
terest. Let us point out that [eg) is doubly-degenerated, the 
hole being located either on the e' g i or on the e' g2 orbitals. 



A. The reference calculation 

The reference calculation includes all effects detailed 
in the preceding section. For the super-conducting com- 
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pound the effective t^g splitting was reported in refer- 
ence^ to be 

AE = e(a lg ) ~ e{e' g ) = 315 meV 

This point corresponds to 9 ~ 61.5° (that is a value of 
9 larger than the one of the regular octahedron 9a ~ 
54.74°) where the crystalline field theory predicts a re- 
verse order between the t 2g orbitals. 

B. Screening effects 

The effect of the screening on the t^g orbital splitting 
can be evaluated by doing a simple CASCI calculation us- 
ing the same fragment, embedding, basis set and orbitals 
as the preceding calculation. Without the screening ef- 
fects, one finds a t2g splitting of 

AE = e(a lg ) - e{e' g ) = 428 meV 

Obviously the screening effects cannot be taken as re- 
sponsible for the qualitative energetic order between the 
aig and e' g orbitals. 

C. Cobalt 3d — oxygen hybridization 

The effect of the hybridization of the cobalt 3d orbitals 
with the neighboring oxygen ligands can be evaluated 
by taking out the oxygen atoms from the quantum clus- 
ter, and treating them as simple —2 point charges at the 
atomic locations. The other parameters of the calcula- 
tion are kept as in the preceding case. The new orbitals 
are optimized at the average-CASSCF level between the 
two |e g ) and the |a lg ) states. It results in a t2 g splitting 
of 

AE = e{a lg ) - s(e' g ) = 40 meV 

for the super-conducting compound. Again the hy- 
bridization of the cobalt 3d orbitals with the neighboring 
oxygens cannot be taken as responsible for the inversion 
of the splitting between the a± g and e' g orbitals. 

D. Long-range electrostatic potential 

The effect of the long-range electrostatic potential can 
be evaluated by restricting the embedding to the NNO 
point charges only, that is to the electrostatic potential 
considered in the crystalline field method. One finds a 
tig splitting of 

AE = e(a lg ) - e(e' g ) = 124 meV 

Once again the results is positive and thus the long-range 
electrostatic potential is not the cause of the crystalline 
field inversion of the ti g splitting. 



E. Orbital radial relaxation 

At this point only few effects on top of the crystalline 
field theory are still treated in the calculation. One of 
them is the radial polarization effect of the 3d orbitals, 
that allows their adaptation to the different occupations 
in the specific |ai g ) and e g ) states. This polarization is 
due to the use of an extended basis set. We thus reduce 
the basis set to a minimal basis set (only one orbital 
degree of freedom per (n, I) occupied or partially occupied 
atomic shell). The minimal basis set was obtained by 
the contraction of the extended one ; the radial part of 
the orbitals being frozen as the one of the the isolated 
Co 4+ high spin state. This choice was done in order to 
keep a basis set as close as possible to the extended one, 
and because only for the isolated atom all 3d orbitals 
are equivalent, and thus have the same radial part. One 
obtains in this minimal basis set a ti g splitting of 

AE = e{a lg ) - e{e' g ) = 41 meV 

At this point we computed the effective orbital energies 
in the sole crystalline field conditions, however the result 
is still reverse than what is usually admitted within this 
approximation. Indeed, the Co 4+ ion was computed in 
the sole electrostatic field of the NNO, treated as —2 
point charges, the calculation is done within a minimal 
basis set, and at the average-CASSCF level. 



F. Further analysis 

In order to understand this puzzling result, we plot- 
ted the whole curve AE{9) (see figure HJ at this level of 
calculation and analyzed separately all energetic terms 
involved in this effective orbital energy difference. 

One sees on figure|2]that the AE(9) curve is not mono- 
tonic, as expected from the crystalline field theory. In- 
deed, while for 9 = the relative order between the a\ g 
and e' g orbitals is in agreement with the crystalline field 
predictions, for 9 = 90° the order is reversed. One should 
also notice that, in addition to the 9q value of the regu- 
lar octahedron, there is another value of 9 for which the 
three ti g orbitals are degenerated. In the physically real- 
istic region of the trigonal distortion (around the regular 
octahedron 9q value) the relative order between the a\ g 
and e' g orbitals is reversed compared to the crystalline 
field predictions. 

Let us now decompose AE{9) into 

• its two-electron part within the 3d shell - 

AE 2 (9) - 

• and the rest referred as 3d single-electron part 

- AE 1 (9). AEi includes the kinetic energy, the 
electron-nucleus and electron-charge interaction, 
and the interaction of the 3d electrons with the in- 
ner shells electrons. 
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FIG. 4: Orbital splitting between the a\ g and e' g orbitals 
when only the nearest neighbor ligands electrostatic field is 
included. The dotted red curve corresponds to the single- 
electron part of the orbital energy difference : AEi, that is 
the kinetic energy (equation {T}), the electron-charge inter- 
action (equation and the interaction with the core elec- 
trons (equation ([3])) . The dashed green curve corresponds to 
the two-electron part of the orbital energy difference : AE2, 
that is the repulsion and exchange terms within the 3d shell 
(equation J!])). The solid vertical line points out the regular 
octahedron 6 value and the dashed vertical line the 9 value 
for the super-conducting compound. 
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where the equations are given in atomic units. Zn refers 
to the nucleus charge of the cobalt atom and the —2 
point charges located at the NNO positions. i?jv is the 
associated electron-charge distance. The sum on \ runs 
over all the orbitals of the cobalt inner-shells. 



Let us now examine the dependence on of each of the 
terms of AE\ and AE 2 ■ 

Kinetic energy : the radial part of each of the 3d or- 
bitals being identical due the the minimal basis set 
restriction, the kinetic part is identical for all 3c? 
orbitals and thus its contribution to AE\ (terms 
labeled [1] of AEi) vanishes. 

Nuclear interaction : obviously this contribution to 
AEi (terms labeled [2] of AEi) strongly depends 
on 6 through the position of the —2 charges. 

Interaction with the inner-shells electrons : this 
term (terms labeled [3] of AEi) depends only on the 
shape of the ti g and inner-shells orbitals. However, 
the minimal basis set does not leave any degree 
of freedom for the relaxation of the inner-shells 
orbital whose shapes are thus independent of 0. 
Similarly, the 3d radial part of the 3d orbitals is 
totally frozen. 

AE2 : finally, the dependence of AE2 can only go 
through the shape of the a\ g and e' g orbitals whose 
radial part is totally frozen due to the use of a min- 
imal basis set. 

If one accepts that the a\ g and e' g orbitals are issued from 
the ti g orbitals of the regular octahedron, their angular 
form is totally given by the symmetry (see eq. E]) and 
both AEi and the third contribution of AE\ should be 
independent of 9. 
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where the x, y and z coordinates are respectively associ- 
ated with the a, b and c crystallographic axes. 

Figure [3] displays both AE\ (dotted red curve) and 
AEi (dashed green curve) contributions to AE. One sees 
immediately that AEi is not at all independent of 9 but 
rather monotonically increasing with 9. It results that 
the above hypotheses of the ti g exclusive origin for the 
e' g orbitals is not valid. Indeed, out of the 9 = 9$ point, 
the only orbital perfectly defined by the symmetry is the 
ai g orbital. The e' g and e g orbitals belong to the same 
irreducible representation (E g ) and can thus mix despite 
the large ti g -e g energy difference. If we name this mixing 
angle a, it comes 



-gi 



e gi cos a + e gi sin a 



-e" ana + e"-cosa 

y L y i 



6 




20 40 60 80 

6° 

FIG. 5: t2 g -e g hybridization angle under the trigonal distor- 
tion. 

Figure [5] displays a as a function of 9. One sees that 
the t2 g -e g hybridization angle a is non null — except 
for the regular octahedron — and a monotonic, increas- 
ing function of 9. Even if very small (±0.6°), this ti g -e g 
hybridization has an important energetic effect, since it 
lowers the the e' g orbital energy while increasing the e g 
one. a is very small but it modulates large energetic fac- 
tors in AE2 : on-site Coulomb repulsions of two electrons 
in the 3d orbitals. The result is a monotonic increasing 
variation of AE2 as a function of 9. The variation of 
the AEi term is dominated by its nuclear interaction 
part and exhibits a monotonic decreasing variation as a 
function of 9, as expected from the crystalline field the- 
ory. The nuclear interaction and t2 g —e g hybridization 
have thus opposite effects on the a\ g -e' g splitting. The 
failure of the crystalline field theory thus comes 
from not considering the t2 g —e g hybridization. 

In the calculations presented in figures Q] and [5l the 
screening effects on the on-site Coulomb repulsions and 
exchange integrals were not taken into account. Thus, 
the absolute value of AE 2 as a function of the hybridiza- 
tion a, is very large and a is very small. When the screen- 
ing effects are properly taken into account, the absolute 
value of AE 2 as a function of a is reduced by a factor 
about 6, and the t2 g -e g hybridization is much larger than 
the values presented in figure [5] Indeed, in the supercon- 
ducting compound, for a realistic calculation including 
all effects, one finds a ~ 13° (6 = 61.5°). 

At this point we would like to compare the a lg -e' g split- 
ting found in the present calculations and the one found 
using DFT methods. Indeed, our splitting (315 meV 
for the superconducting compound) is larger than the 
DFT evaluations (always smaller < 150 meV). This 
point can be easily understood using the single-electron 
and two-electron part analysis presented above. In- 
deed, while the single-electron part is perfectly treated 
in DFT calculations, the two-electron part is treated 
within the exchange-correlation kernel. However these 
kernels are well known to fail to properly reproduce the 



strong correlation effects present in the transition metal 
opened 3d shells. One thus expect that while the single- 
electron part of the atomic orbital energies is well treated, 
the two-electron part is underestimated, resulting in an 
under-evaluation of the a\ g -e' g splitting, as can be clearly 
seen from figure [5] 



IV. OTHER CASES 

We considered up to now a Co 4+ ion, that is five elec- 
trons in the 3d shell, and a fixed metal-ligand distance, 
-Rm-o- Let us now examine the effect of the distance 
KM — O and the band filling on the a\ g -e' g splitting. The 
calculations presented in this section follow the same pro- 
cedure as in sections IIIIE1 IIIIFI For different fillings a 
typical example in the transition metal oxides family was 
used to define the type of metallic atom and metal oxygen 
distances. Minimal basis set issued from full contraction 
of the basis set given in reference^ will be used. 



A. The effect of the Co— O distance 

Figure [5] displays the a\ g -e' g energy splitting as a func- 
tion of the distortion angle 9 and for different distances. 
The range of variation : from 1.8A to 1.95A, includes 
all physically observed distances in C0O2 layers. One 
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FIG. 6: Orbital splitting between the a\ g and e' g orbitals 
for a 3d 5 transition metal and for different metal-ligand dis- 
tances. Only the nearest neighbor ligands electrostatic field is 
included in the calculation. The dotted red curve corresponds 
to Rco-o = 1.8 A, the solid black curve corresponds to the 
superconducting compound (Rco-o = 1.855 A), the magenta 
dashed curve corresponds to Rco-o = 1-9 A, and finally the 
dot-dashed blue curve corresponds to Rco-o = 1-95 A. 

sees immediately that despite the large variation of the 
metal-ligand distance, the relative order of the a\ g and 
e' g orbitals remains identical. The main effect of KM — O 
is thus to renormalize the amplitude of the splitting, low- 
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ering the splitting for larger distances and increasing it compared to the crystalline field predictions. As for the 
for smaller ones. 3d 5 system, there is a slight G g -e g hybridization that is 



B. 3d 1 



responsible for the i 2g orbitals order. 



The simplest filling case corresponds to only one elec- 
tron in the 3d shell. This is, for instance, the case 
of the NaTiC>2 compound. The calculations were done 
using the average Ti-0 distance found in NaTiO^ : 
flxi-o = 2.0749A. 

In this case, AE 2 = and AE(9) = AE 1 (9) behaves 
as pictured in figured] The a\ g orbital is of lower energy 
than the e' g for 9 > 9o and of higher energy for 9 < 6$. 
This result is in perfect agreement with the crystalline 
field theory. 
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A simple example of the 3d 2 filling in transition metal 
oxides is the LiV02 compound. Indeed, the vanadium 
atom is in the V 3+ ionization state. We thus used a 
metal oxygen distance of Ry-o = 1.9787A^. Figure 
displays the a\ g -e' g splitting as well as its decomposition 
into the single-electron and two-electron parts. As in the 
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D. 3d J 

Examples of 3d 3 transition metal oxides are found eas- 
ily in the chromium compounds. Let us take for instance 
the NaCrC>2 systemM. The metal oxygen distance is 
thus : i?o-o — 1.9Q1A. Figure [8] displays the ai g -e' g 
orbital splitting as well as its decomposition into single- 
and two-electron parts. As usual the single-electron part 
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FIG. 8: Orbital splitting between the a\ g and e' g orbitals for a 
3d 3 transition metal. Only the nearest neighbor ligands elec- 
trostatic field is included in the calculation. The dotted red 
curve corresponds to the single-electron part of the orbital 
energy difference : AEi, that is the kinetic energy (equa- 
tion {TJ), the electron-charge interaction (equation (J2J) and 
the interaction with the core electrons (equation ©) . The 
dashed green curve corresponds to the two-electron part of 
the orbital energy difference : AE2, that is the repulsion and 
exchange terms within the 3d shell (equation (0J)). 



FIG. 7: Orbital splitting between the ai g and e' g orbitals for a 
3d 2 transition metal. Only the nearest neighbor ligands elec- 
trostatic field is included in the calculation. The dotted red 
curve corresponds to the single-electron part of the orbital 
energy difference : AEi, that is the kinetic energy (equa- 
tion {TJ), the electron-charge interaction (equation ©) and 
the interaction with the core electrons (equation (O) . The 
dashed green curve corresponds to the two-electron part of 
the orbital energy difference : AE2, that is the repulsion and 
exchange terms within the 3d shell (equation Q). 

3d 5 case (figure 2]), the single-electron and two-electron 
parts behave in a monotonic way as a function of 9, and 
in an opposite manner. In the present case, however, the 
two-electron part always dominates over the one-electron 
part and the a\ g -e' g orbital splitting is always reversed 



and the two-electron part are monotonic as a function of 
6 but with slopes of opposite signs. This case is quite 
similar to the 3d 5 case since none of the single- and two- 
electron parts dominates the t 2g orbital splitting over the 
whole range. Indeed, for small values of 9, the crystalline 
field effect dominates and the a\ g orbital is above the e' g 
ones while, for large values of 9, the two-electron part 
dominates and the a± g orbital is again above the e' g ones. 
In a small intermediate region the order is reversed. In 
the realistic range of 9 (9 ~ 9q) there is a strong competi- 
tion between the two effects (quasi-degeneracy of the a\ g 
and e' g orbitals) and no simple theoretical prediction can 
be made. The crystalline field theory is not predictive 
but the present calculations cannot be considered as pre- 
dictive either, since all the neglected effects may reverse 
the aig-e' order. 
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V. DISCUSSION AND CONCLUSION 

In the present work we studied the validity of the crys- 
talline field theory under the application of a trigonal 
distortion on the regular octahedron. Under such a dis- 
tortion, the T 2g irreducible representation (irrep) of the 
Oh group spits into A\ g and E g irreps (T 2g — > A\ g ®E g ), 
while the e g irrep remains untouched [E g — ► E g ). The 
hybridization between the t 2g and e g orbitals thus be- 
come symmetry allowed, even if hindered by energetic 
factors. This hybridization is not taken into account in 
the crystalline field theory. It is however of crucial im- 
portance for the relative order between the former t 2g 
orbitals and the reason of the failure of the crystalline 
field theory to be predictive. Indeed, due to the t 2g -e g 
orbitals hybridization, the two-electron part of the e' g or- 
bital energy becomes dependant of the amplitude of the 
distortion and of opposite effect to the single-electron 
part. The relative order of the t 2g orbitals thus depends 
on the competition between these two effects and as a 
consequence of the band filling. 

In this work we studied the Oh to D^d distortion, how- 
ever one can expect similar effects to take place for other 
distortions of the regular octahedron. The condition for 
these effects to take place is that the T 2g irreducible rep- 
resentation splits into a one-dimensional irrep (A) and 
the same two-dimensional irrep (E) as the one the e g 
orbitals are transformed to 

T 2g — > A®E 

Eg * E 

Indeed, under such a distortion, t 2g -e g hybridization phe- 
nomena are allowed. The distortion should thus trans- 
form Oh into sub-groups that keep the C3 (111) sym- 
metry axis : C3, C3,,, D3, Sq and D^d- Examples of 



such deformations are the elongation of the metal-ligand 
distance of one of the sets of three symmetry related lig- 
ands, or the rotation of such a set three ligands around 
the (111) symmetry axis. For instance, one will expect 
that t 2g -e g hybridization will also take place in trigonal 
prismatic coordination. 

However, in real systems like the sodium cobaltites, 
these distortion do not usually appear alone but rather 
coupled. For instance, in the squeezing of the metal layer 
between the two oxygen layers observed as a function of 
the sodium content in Na x Co02, the Co-0 bond length 
and the three-fold trigonal distortion are coupled. Since 
this composed distortion belongs to the above-cited class, 
the t 2g -e g hybridization will take place and the relative 
orbital order between the a\ g and e' g orbitals will be qual- 
itatively the same as in figure [4j The bond length mod- 
ification at equal distortion angle, 9, will only change 
the quantitative value of the orbital splitting, but not 
its sign. A bond elongation reduces the splitting a bond 
compression increases it. One can thus expect in sodium 
cobaltites that the a\ g -e' g orbital energy splitting will de- 
crease with increasing sodium content. The reader should 
however have in mind that the effects of this splitting 
reduction will remain relatively small compared to the 
band width as clearly seen in reference^. In fact, one 
can expect that a large effect will be the modification 
of the band dispersion due not only to the bond length 
modification, but also to the t 2g -e g hybridization. 
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